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ABSTRACT 

A sample of 20 powerful extended radio galaxies with redshifts between zero 
and two were used to determine constraints on global cosmological parameters. 
Data for six radio sources were obtained from the VLA archive, analyzed, and 
combined with the sample of 14 radio galaxies used previously by Guerra & Daly 
to determine cosmological parameters. The new results are consistent with our 
previous results, and indicate that the current value of the mean mass density 
of the universe is significantly less than the critical value. A universe with Vt m 
in matter of unity is ruled out at 99.0% confidence, and the best fitting values 
of Vt m in matter are O.lOloiio an d — 0.25t[|;25 assuming zero space curvature and 
zero cosmological constant, respectively. 

The radio properties of each source are also used to determine the density of 
the gas in the vicinity of the source, and the beam power of the source. The six 
new radio sources have physical characteristics similar to those found for the 
original 14 sources. The density of the gas around these radio sources is typical 
of gas in present day clusters of galaxies. The beam powers are typically about 
10 45 erg s _1 . 
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1. Introduction 

The future and ultimate fate of the universe can be predicted given a knowledge of 
the recent expansion history of the universe (assuming the universe is homogeneous and 
isotropic on scales greater than the current horizon size). This recent expansion history 
can be probed by studying the coordinate distance to sources at redshifts of one or two; 
the coordinate distance is equivalent to the luminosity distance or angular size distance 
multiplied by factors of (1 + z). The advantage of determining cosmological parameters 
using the coordinate distance (luminosity distance or angular size distance) is that this 
distance depends on global, or average, cosmological parameters. It is independent of the 
way the matter is distributed spatially, of the power spectrum of density fluctuations, of 
whether the matter is biased relative to the light, and of the form or nature of the dark 
matter (assuming that the universe is homogeneous and isotropic on large scales). 

It was shown in 1994 that powerful double-lobed radio galaxies provide a modified 
standard yardstick that can be used to determine global cosmological parameters (Daly 
1994, 1995), much like supernovae can be used as modified standard candles. The method 
was applied and discussed in detail by Guerra & Daly (1996, 1998), Guerra (1997), and 
Daly, Guerra, & Wan (1998) who found that the data strongly favor a low density universe; 
a universe with Q m = 1 was ruled out at 97.5 % confidence. 

It was shown by Daly (1994, 1995) that the radio properties of these sources could be 
used not only to study global cosmological parameters, but also to determine the ambient 
gas density, beam power, Mach number of lobe advance, and ambient gas temperature of the 
sources and their environments. The characteristics of the sources and their environments 
are presented and discussed in a series of papers (Wellman & Daly 1996a,b; Wan, Daly, & 
Wellman 1996; Daly 1996; Wellman, Daly, & Wan 1997a,b; Wan & Daly 1998a,b; Wan, 
Daly, & Guerra 1998). 

Radio maps of six powerful double-lobed radio galaxies were extracted from the Very 
Large Array (VLA) archives at the National Radio Astronomy Observatory (NRAO), and 
analyzed in detail. New results on global cosmological parameters, ambient gas densities, 
and beam powers are presented here. 

The expanded sample is described in §2. The new results are presented in §3. The 
implications of the results are discussed in §4. 
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2. Expanded Sample 

Each powerful extended radio galaxy (also known as a "classical double") has a 
characteristic size, -D*, that predicts the lobe- lobe separation at the end of its lifetime (Daly 
1994, 1995; Guerra & Daly 1998). The parameters needed to compute the characteristic 
size are the lobe propagation velocity, vl, the lobe width, ol, and the lobe magnetic field 
strength, B^. These three parameters can be determined using radio maps with arc-second 
resolution at multiple frequencies, such as those produced with the VLA or MERLIN 
(Multi-Element Radio Linked Interferometer Network). Multiple-frequency data are needed 
to use the theory of spectral aging to estimate the lobe propagation velocity (e.g., Myers & 
Spangler 1985). In addition, these maps must have the necessary angular resolution and 
dynamic range to image sufficient portions of the radio bridges. 

Two published data sets, Leahy, Muxlow, & Stephens (1989) and Liu, Pooley, & Riley 
(1992), have radio maps of powerful extended radio galaxies at multiple frequencies which 
are sufficient to compute all three parameters used in determining D*. These data were 
used to compute for 14 radio galaxies (Guerra & Daly 1996, 1998; Guerra 1997; Daly, 
Guerra, & Wan 1998). Current efforts to expand the data set are underway, and include 
searches through the VLA archive. The VLA archive search has already yielded the desired 
data for six radio galaxies, and new results including these six sources are presented here. 

Data were selected from observations of powerful extended radio galaxies from the 3CR 
sample (Bennett 1962) on the basis of the observation frequency and array configuration 
used. An observation in the VLA archive was a candidate if the lobe-lobe angular size of 
the source was 10 to 40 times the implied beam size, and four hours separated the first 
and last scans. These observations should resolve the source sufficiently and have enough 
Mf-coverage to image the bridges. From candidate observations at both L and C band, we 
have successfully imaged data from six sources. The VLA archive data sets used here are 
listed in Table [TJ. 

Radio imaging was performed using the NRAO AIVS software package. The uv data 
needed minimal editing, and initial calibration was performed in the standard manner using 
AIVS. As diagnostics, initial images were made using the AIVS task IMAGR both without 
and with the CLEAN deconvolution algorithm. The final radio maps were produced with the 
SCMAP task, which performs self-calibration along with the IMAGR and CLEAN tasks. 

Parameters used in AIVS for imaging, deconvolution, and calibration are chosen to 
produce radio maps for a given source that can easily be used for spectral aging analysis. 
This is particularly important where observations at different frequencies are produced by 
different observers. Although some of these radio maps exist in published form, different 
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choices of parameters (such as the restoring beam) are often used to produce these radio 
maps. Thus, reducing the raw uv data insures that the data analysis can be performed 
consistently between radio maps. 

The computation of a^, Bl, and vi were performed here in the same manner as 
Wellman, Daly, & Wan (1997a,b). The deconvolved lobe width, a^, is measured 10/t -1 
kpc from the hot spot toward the host galaxy. The lobe magnetic field strength, B^, is 
computed from the deconvolved surface brightness and bridge width measured lO/i" 1 kpc 
from the hot spot toward the host galaxy. The lobe propagation velocity, vl, is computed 
on the basis of spectral aging along the imaged portions of the bridge, and the magnetic 
field used in spectral aging is computed using values measured 10/i _1 kpc and 2hh~ l kpc 
from the hotspot. In this paper, all parameters are computed with b = 0.25 which means 
magnetic fields are computed to be 0.25 times the minimum energy values, and without an 
a — z correction which refers to a correction related to the observed correlation between 
spectral index and redshift (see Wellman 1997; Wellman, Daly, Wan 1997a,b; Guerra 1998; 
Guerra & Daly 1998 for details). It was found by Guerra & Daly (1998) that constraints 
on cosmological parameters did not depend on these choices, and very similar results are 
obtained independent of the value of b and of whether an a — z correction is applied. 

These three parameters are computed for each bridge, with the exception of one bridge 
in 3C 324 which was not imaged along its length sufficiently The characteristic core-lobe 
size, r*, was computed for each bridge using using the equation (see Guerra & Daly 1998): 



2/3/3 



r*oc I-— J v\~ m . (1) 
\B L a L J 

The characteristic (lobe- lobe) size, is taken to be the sum of both r* (or in the case 
of 3C 324, D* = 2r#), and is normalized so that of Cygnus A (3C 405), a very low 
redshift source in our sample, is equal to the average size of the full population of powerful 
classical double radio galaxies at very low redshift. Constraints on cosmological parameters 
determined using this method are independent of the normalization of D*. As discussed in 



§ |3.1| , the best fit value of the one model parameter (3 is determined simultaneously with the 
best fit values for the two cosmological parameters that enter, Q m and Q^, the normalized 
values of the current values of the mean mass density and the cosmological constant (see 
Guerra & Daly 1998; Daly, Guerra, & Wan 1998). 

Table [| presents the six new values in the last column, assuming the best fit value 



of (5 = 1.75 (see § |3.1|) . Source name and redshift are listed in the first two columns. The 
third column lists the redshift bin corresponding to the assignments in Guerra & Daly 
(1998) and Table |3| below. The lobe-lobe angular size of the source is listed in the fourth 
column, and the fifth and sixth columns list the core- lobe characteristic sizes, r*. 
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3. Updated Results 

3.1. Constraints on Cosmological Parameters 

The subsample of sources with estimates of the characteristic size, D*, has been 
increased from 14 to 20, as discussed above in §^|. Only sources with physical sizes, defined 
as the projected separation between the radio hot spots, greater than 20 h^ 1 kpc can 
be used to determine a characteristic size. This is because smaller sources are typically 
not sufficiently resolved that the radio data are useful, and the radio lobes of smaller 
sources are interacting with the interstellar medium of the host galaxy rather than the 
intergalactic/intracluster medium. It was decided that this same criterion should be applied 
to the larger comparison sample of powerful 3CR radio galaxies. Thus, the sample of radio 
galaxies used to determine the redshift evolution of the physical size has been reduced from 
82 to 70; twelve radio galaxies were cut because the physical separation between their lobes 
was less than 20 h~ x kpc. This has a rather small impact on the actual means and standard 
deviations of the parent population, as can be seen by comparing Table ^| of this paper 
with Table 1 of Guerra & Daly (1998). The average lobe-lobe separations as a function 
of redshift are listed in Table |3| for three example choices of cosmological parameters 
(matter-dominated, curvature-dominated, and spatially flat with non-zero cosmological 
constant). 

To solve simultaneously for the model parameter (3 and the cosmological parameters 
Q m and Q\, the ratio of for each source to (D), the average lobe-lobe size of the 
parent population in the corresponding redshift bin, is fit to a constant, independent of 
redshift. The value of the constant is a free parameter, so the normalization of Cygnus 
A does not affect the fits in any way. Figure [1] illustrates the cosmological dependence of 
(D)/D* on the coordinate distance (a a r). For (3 = 1.75, (D)/D* is proportional to (a r) 16 ; 
thus ((D) I D*)(a r)~ ls is independent of cosmological parameters. The data can be 
compared to several different sets of cosmological parameters on a single figure by plotting 
((D) / D*)(a r)~ ls for each data point and comparing this with (a r) -1 ' 6 curves obtained 
for different sets of cosmological parameters, as is shown in Figure [l]. 

The hypothesis is that, for the correct choice of cosmological parameters, (D) / D* = 1, 
so that the values of ((D) / D*)(a r)~ ls for all 20 radio galaxies should follow a curve that, 
at each z, is parallel to the curve (a r)~ 1M obtained for that particular choice of cosmological 
parameters. Figure [I], shows ((D) / D 1r )(a Q r)~ l & for the the six new points and 14 original 
points as a function of z (the six new points are denoted with stars). Also drawn on this 
figure are curves of (a D r)~ L6 obtained for specified values of cosmological parameters. In 
this figure, all of the curves pass through Cygnus A, though this is not required when we 
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actually solve for best fitting cosmological parameters. It is clear that curves obtained for a 
low density universe describe the data points quite well, and the curve describing a universe 
with Q m = 1 does not follow the data points. Note that actual fits are performed with a 
free normalization that does not require the predicted curves to pass through Cygnus A, 
the lowest redshift point. Figure [l] is for illustrative purposes and is not part of any fits. 

For all 20 sources, the chi-squared for fitting the ratio to a constant is computed for 
relevant values of (3, Q m , and Q\. It is found that the best fit value of /3 is (3 — 1.75 ± 0.25. 
This result is insensitive to the choice of Q m and Q\, and there appears to be no significant 
covariance between (3 and cosmological parameters (see Figures 0a,b). 

The confidence contours in the Q m - plane are shown in Figures [3] and f|. The 
probability associated with a given range of Q m and fl\ independent of j3 is shown in Figure 
|3] (referred to as two-dimensional confidence intervals). In Figure |], the projection of a 
confidence interval onto either axis (fl m or Q\) indicates the probability associated with 
a given range of that one parameter, independent of all other parameter choices (referred 
to as one-dimensional confidence intervals). Both figures illustrate how this method and 
the data are most consistent with a low density universe; Q m ^ 0.15 with 68% confidence, 
^ 0.5 with 90% confidence, and Q m < 1.0 with 99% confidence. The constraints on Qa 
are not as strong, and values of Q\ from zero to unity are consistent with the data. 

The best fit value of (3 = 1.75 ± 0.25 is consistent with the previous estimates of 
(3 ~ 1.5 ± 0.5 (Daly 1994), and (3 ~ 2.1 ± 0.6 (Guerra & Daly 1996, 1998), but with 
significantly reduced uncertainties. Similarly, the constraints on cosmological parameters 
are consistent with previous estimates, (Daly 1994; Guerra & Daly 1996, 1998; Guerra 1997; 
Daly, Guerra, & Wan 1998) but with smaller error bars. It is apparent in Figures ^] and |] 
that these data and method strongly favor a low density universe; a universe where fl m = 1 
is ruled out with 99.0 % confidence independent of and (3 . This will be discussed in 
more detail in §|| 



3.2. Ambient Gas Densities and Beam Powers 

Daly (1994, 1995) showed that the radio properties of a powerful extended radio source 
could be used to estimate the beam power, Lj, and density of the gas in the vicinity of the 
source, n a . The method is described in detail by Wellman (1997), Wellman, Daly, & Wan 
(1997a), Wan (1998), and Wan, Daly, & Guerra (1998). Values for the ambient gas density 
and beam power for the 14 sources in the original sample are described in these papers. 
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The basic equations are: 

Lj oc a\ B 2 v L , (2) 

and 

B 2 

n a oc — f- ; (3) 

the normalizations are given in the references listed above. The values obtained for the six 
new radio galaxies in our sample are listed in Table [|. Also listed in Table f| are all the 
input parameters used to compute Lj, n a , and D*. 



4. Discussion 

A parent population of 70 powerful extended classical double radio galaxies with 
redshifts between zero and two was used to define the evolution of the mean or characteristic 
size of these sources as a function of redshift. An independent estimate of the mean or 
characteristic size of a given source was possible for a subset of 20 of these radio galaxies 
for which extensive multiple frequency radio data was available. Requiring that the two 
measures of the mean source size have the same redshift behavior allows a simultaneous 
determination of three parameters: the one model parameter (3, and the two cosmological 
parameters Vt m and Q\ (assuming that the only significant cosmological parameters today 
are the mean mass density, a cosmological constant, and space curvature). 

The method was applied to this data set, and interesting new constraints are presented. 
It is found that the model parameter is very tightly constrained to be (3 = 1.75 ± 0.25 (see 
Figure H), consistent with previous estimates, and that this model parameter is independent 
of cosmological parameters. For a value of (3 = 1.75, the characteristic source size is 

oc [B^a^ 1 ^ vj}l 12 , which indicates that it is necessary to have multiple-frequency 
radio data in order to estimate D*, owing to its vl dependence. 

The data strongly favor a low density universe; a universe with Q m = 1 is ruled out 
with 99% confidence, independent of the value of Q\ or (3. Either space curvature or a 
cosmological constant, or both, are allowed. The main conclusion is that Q m is low, but, 
at this point, the method and data do not allow a discrimination between whether space 
curvature or a cosmological constant is important at the present epoch. 

It is interesting to note that the lowest reduced chi-squared obtained is 0.96 for 
Q m ~ —0.25 and fl^ ~ 0. This value is slightly greater and closer to unity than the 
minimum reduced chi-squared obtained by Guerra & Daly (1998), which indicates that this 
sample of 20 sources has a reasonable distribution around any model predictions. This 
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convergence to unity with increasing sample size suggests that this method and its statistics 
are reliable. 

The best fit for cosmological parameters in the physically relevant half-plane of fl m > 
is fl m = and Q\ = 0.45 with a reduced chi-squared of 0.98 (also close to unity). However, 
Figures § and ^ clearly show that our results are still consistent Q\ = 0. 

The radio data also allow a determination of the density of the ambient gas in the 
vicinity of each radio source, and the beam power of each source; the values of these 
quantities are presented. The sources lie in high-density gaseous environments like those 
found in low-redshift clusters of galaxies. Typical beam powers for the sources are 
~ 10 45 erg s -1 . 

Special thanks go to Rick Perley for his aid in extracting data from the VLA archive. 
The authors would also like to thank Katherine Blundell, Chris Carilli, Miller Goss, Paddy 
Leahy, Wil van Breugel, and Dave Wilkinson for helpful discussions. It is a pleasure to 
acknowledge Ingrid Stairs, Rick Balsano, and other members of the pulsar group who 
graciously allowed us to use their computer facilities. This work was supported in part by 
the U.S. National Science Foundation, the Independent College Fund of New Jersey, and a 
grant from W. M. Wheeler III. 
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Table 1. VLA Archive Data Sets 



Source 


Band 


Program ID 


Listed 
Observer 


Frequency 
(GHz) 


Array 
Config. 


Date 


3C 244.1 


L 


POOL 


Pooley, G. 


1.411 


B 


09/17/82 




C 


AF213 


Fernini, I. 


4.872 


B 


12/23/91 


3C 337 


L 


AR123 


Rudnick, L. 


1.452 


A 


02/21/85 




C 


AR123 


Rudnick, L. 


4.885 


B 


06/01/85 


3C 325 


L 


AV153 


van Breugel, W. 


1.465 


A 


12/05/88 




C 


AF213 


Fernini, I. 


4.885 


B 


12/23/91 


3C 194 


L 


AV164 


van Breugel, W. 


1.465 


A 


05/11/90 




C 


AV164 


van Breugel, W. 


4.885 


A 


05/11/90 


3C 324 


L 


AR123 


Rudnick, L. 


1.452 


A 


02/21/85 




C 


AR123 


Rudnick, L. 


4.885 


B 


06/01/85 


3C 437 


L 


AV164 


van Breugel, W. 


1.465 


A 


05/11/90 




C 


AV164 


van Breugel, W. 


4.885 


B 


04/21/89 



Table 2. Radio Galaxies with D± Presented Here 



Source 


z 


Bin 


e 

(arcsec) 


(h- 1 kpc) 


' * 

(h- 1 kpc) 


(h- 1 kpc) 


3C 244.1 


0.43 


2 


50.8 


151 ± 16 


135 ± 16 


286 ± 22 


3C 337 


0.63 


3 


43.5 


150 ± 24 


65 ±7 


214 ± 25 


3C 325 


0.86 


3 


15.8 


164 ± 52 


66 ±14 


230 ± 54 


3C 194 


1.19 


4 


14.2 


105 ± 22 


92 ± 15 


197 ± 27 


3C 324 


1.21 


5 


10.2 


75 ± 18 




149 ± 35 


3C 437 


1.48 


5 


36.7 


55 ±6 


48 ±5 


103 ±8 


a Computed 


assuming 


; P = L75, fi 


= 0.1, ft A 


= 0, b = 0.25 and not including a 


— z correction. 



'r* for only one bridge. 



Table 3. The Average Lobe-Lobe Sizes for Powerful 3CR Radio Galaxies. 



(D) (h- 1 kpc) 



Bin 


z Range 


Sources 


n a = 1.0, Q A = 0.0 


n o = 0.1, tt A = 0.0 


n Q = 0.1, ft A = 0.9 


1 


0.0-0.3 


3 


66±14 


68±14 


72±13 


2 


0.3-0.6 


13 


202±45 


224±50 


259±57 


3 


0.6-0.9 


23 


148±17 


174±20 


209±24 


4 


0.9-1.2 


16 


107±24 


133±29 


165±36 


5 


1.2-1.6 


9 


91±32 


122±43 


152±53 


6 


1.6-2.0 


6 


66±19 


92±26 


114±33 
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Table 4. Summary of New Source Properties 



Source 


z 


a L a 


-Bio b 


-E>25 C 


v L d 


n a e 


logV 


3C 244.1 


0.43 


5.5 ±0.3 


3.5 ±0.2 


2.5 ±0.2 


1.1 ±0.2 


1.0 ±0.4 


44.35 ±0.10 






4.7 ±0.4 


5.3 ±0.3 


3.2 ±0.2 


1.7 ±0.3 


1.0 ±0.4 


44.75 ±0.11 


3C 337 


0.63 


6.2 ±0.5 


3.5 ±0.2 


2.3 ±0.1 


1.5 ±0.5 


0.6 ±0.4 


44.57 ±0.15 






10.5 ±0.4 


3.8 ±0.2 


2.8 ±0.2 


1.1 ±0.2 


1.2 ±0.4 


44.99 ±0.10 


3C 325 


0.86 


3.3 ± 1.1 


6.8 ±0.9 


4.0 ±0.3 


2.0 ±0.5 


1.2 ±0.6 


44.73 ± 0.23 






4.4 ±0.8 


12.7 ± 1.0 


4.6 ±0.3 


3.0 ±0.8 


1.9 ± 1.0 


45.70 ±0.17 


3C 194 


1.19 


4.5 ± 1.0 


7.1 ±0.6 


5.0 ±0.3 


1.8 ±0.3 


1.6 ±0.6 


45.00 ±0.16 






5.6 ±0.8 


6.7 ±0.5 


5.6 ±0.4 


2.0 ±0.3 


1.2 ±0.4 


45.16 ±0.13 


3C 324 


1.21 


4.5 ±0.9 


10.3 ±0.9 




2.1 ±0.7 


2.4 ±1.5 


45.38 ±0.19 


3C 437 


1.48 


11.4 ±0.8 


6.9 ±0.4 


4.2 ±0.2 


4.6 ±0.8 


0.2 ±0.1 


46.18 ±0.10 






10.5 ±0.8 


9.4 ±0.6 


6.7 ±0.4 


6.0 ±1.0 


0.2 ±0.1 


46.49 ±0.10 



a Lobe radius, 10 h kpc behind hot spot, in h kpc. 

b Minimum energy magnetic field, 10 h^ 1 kpc behind hot spot, in /i 2 / 7 10~ 5 G. 
c Minimum energy magnetic field, 25 h^ 1 kpc behind hot spot, in /i 2 / 7 10~ 5 G. 
d Lobe advance speed, in 10~ 2 c. 
e Ambient gas density, in 10 _3 /i 1 / 2 cm -3 . 

f Logarithm of the luminosity in directed kinetic energy, in /i~ 2 erg s _1 . 
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Fig. I.— The quantity ((D)/£>,)(a r) -:L6 , computed assuming (5 = 1.75, b = 0.25 and not 
including a — z correction. This measured quantity is independent Q m and Q\. Open Stars 
denote the six new data points from the VLA archive. For different choices of Q m and ^a, 
the predicted redshift evolution of (a D r)~ L6 is plotted for comparison. 



- 15 - 




Fig. 2. — The 68% and 90% confidence intervals for (5 and Q m , where (a) f2A=0 (no 
cosmological constant)... 




Fig. 2.— ...and (b) fi A = 1 



— Q m (spatially flat universe). 
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Fig. 3. — The 68% and 90% confidence intervals for ranges of both Q m and f2 A , independent 
of (3. (Two-dimensional) 
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Fig. 4. — The projections of the 68%, 95%, and 99% confidence intervals onto either axis (Q m 
or Q\) indicates the probability associated with the range in that one parameter, independent 
of all other parameter choices. (One-dimensional) 



